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Antenna Coupling Error in Direction Finders * 

Charles W. Harrison, Jr. 

(August 15, 1960; revised January 19, 1961) 

Antenna coupling or scattering error in an interferometer-angle-measuring system 
consisting of two identical base-loaded quarter-wave vertical antennas over a perfectly 
conducting earth is investigated. A curve is supplied for the error in angle as a function of 
antenna spacing for an incident vertically polarized electric field arriving in the plane of 
the antennas, when the ratio of antenna length to radius is 75 and the base loads are resistors 
of 50 ohms. 



1. Introduction 

An interferonieter-angle-measuring system for 
missile trajectory determination — in rudimentary 
form — consists of four identical base-loaded vertical 
antennas, one-quarter wave or less in length, creeled 
at the vertices of a square. The phase angles be- 
tween the voltages developed across identical load 
impedances in diagonal antenna pairs are measured. 
These angles establish the direction cosines of the 
missile from which the signal emanates. In addition 
to electronic missile tracking equipment, many fixed 
antenna direction finders utilize this principle. 

There are several sources of error in interferometer- 
angle-measuring systems. It is convenient to class 
these errors under two headings: geometrical and 
electrical. Simmons [1959] has studied geometrical 
errors in a radio interferometer. Under electrical 
errors, Carswell and Flammer [1957] investigated the 
error clue to multipath. In the present paper antenna 
coupling error in a two-antenna interferometer is 
investigated theoretically and numerically. The 
model employed consists of two isolated parallel 
nonstaggered antennas of the same length and 
radius, center-loaded by identical impedances. The 
incident electric field is assumed to be polarized 
parallel to the axes of the antennas, but this restric- 
tion in the theory may be removed. The azimuth 
angle of wave arrival that is indicated by the equip- 
ment is computed and compared to the assumed 
angle of wave arrival. The computed error angle 
then applies to a base-loaded two-antenna array 
over a perfectly conducting earth for which the 
loading impedances are halved. Second order 
accuracy is achieved in the theoretical development 
and in the numerical results of this report. In 
another paper Harrison [1960] studied the more 
general case of the four-antenna direction finder and 
obtained results accurate to zeroth order. 

Some insight into the nature of antenna coupling 
error may be gained from the following remarks: 
Each antenna in the array continuously samples the 
scattered field from the other antennas comprising 
the array. The resultant field acting on each antenna 
in the interferometer is the vector sum of the field 
maintained along its surface by the distant telemetry 
transmitter on the missile and the several scattered 



fields maintained by the currents in the other 
antennas. Thus the coordinates of the missile, 
determined by the phase-comparison method, must 
be in error. The error is a, function of the polarization 
of the electric field, the coordinates of the signal 
source, the antenna dimensions (length and radius), 
the spacing between elements, the load impedance, 
and the properties of the earth. In this paper the 
scattering error in an interferometer that is far from 
the earth, or that is erected on a very highly conduct- 
ing plane, is determined. The effect of an imper- 
fectly conducting earth is not treated. 

2. Integral Equations for the Currents in an 
Unloaded Two-Element Receiving Array 
and Their Solutions 



"""Consider two identical non-staggered parallel 
receiving antennas, labeled 1 and 2, oriented as shown 
in figure 1 with respect to a system of rectangular 



1 Contribution from the Sandia Laboratory, Albuqerque, N. Mex. 




Figure 1. Two-element receiving antenna array (load terminals 
short-circuited) . 
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and spherical coordinates with common origin. The C h Tn/ , XT ^ n( /N 7 , . 47r r ^„ A , TTrt1 

half-length of each antenna is A, the radius a, and J_^ 7 ^ > if ^ z ) dz =~3 J ^ C ™s(3z+U a } 
the spacing between them is b. It is assumed that 
A>>a, and /3a<<Cl. Here ft is the radian wave 
number. The terminals of the receiving elements 

are located at £=0, y= ±— , and 2=0. The incident 

z 

electric field E arrives from the distant source at 
the azimuth angle 3>, measured from the positive 
x-axis, and is directed parallel to the axes of the 
wires. The currents in antennas 1 and 2, when the 
terminals are short-circuited, are designated I\ E {z) 
and I 2E (z), respectively, and are functions of E, <£, 
and the electrical dimensions of the array /3a, ph, 
and ftb. From symmetry considerations it is clear 
that I\e{z)9 £ I 2E {z) except when $=mr, n=0, 1, 
2 . . . . 

When the electric field is directed parallel to the 
axes of the wires, the simultaneous integral equations 
for determining the currents I\ E {z) and I 2E {z) are 
[King, 1956, p. 517, eqs (2a) and (2b)] 



f* I 1E (z')K a (z,z>)dz' + f* I 2E {z')K,{z,z')dz' 

J -h J -k 

= -i^/aco S /3 2 +LV? s,n *|, 

J -h J -h 

= -Jy{acos / 3z+^- / >*}. 



(1) 



Here 



e -& V(~ 



-z') 2 +a 2 



V(2-s') 2 + « 2 



(2) 



(3) 



K b {z,z') is obtained from (3) by writing b for a, 
f is the characteristic resistance of space, C x and 
C 2 are constants, and 

U=-E/$. (4) 
The currents I 1E (z) and I 2E (z) may be written 

Ii E (z)=P E (z)+mz) (5) 

I 2E (z)=I s E (z)-I E (z) (6) 

where I%(z) and I E (z) are the symmetrical and 
antisymmetrical components of current in either 
antenna. The assumed directions of these currents 
are indicated in figure 1. 

Substituting (5) and (6) into (1) and (2), the 
following independent integral equations for the 
currents I%(z) and I E {z) are obtained: 



/. 



* I%(z')K°(,z,z')dz' = -j~ {C s cosj3z+U°} (7) 
n J 



In (7) and (8), 

K'(z,z')=K a (z,z')+K b (z,z'), 
K a (z,z')=K a (z,z')-K b (z,z'), 

U S =U cos (^ sin *\ 
0b 



U a =jU sin (^sin $) 






&= 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 

(14) 



The integral equation for the unloaded current in 
an isolated receiving antenna consisting of a single 
conductor is identical in form to (7) and (8). It 
follows that solutions for the currents I E (6) and 
I E (o) may be obtained from investigations of the 
receiving dipole by appropriate substitutions for 
the kernel and the [/-function. These currents 
may also be obtained from the relations 



/i(0) = - 



20h s e U s 



2ph s e U cos 0j sin $) 



Zsl\£l 2 



and 



I%(0) = 



2ph a e U a 



2$h a e U sin Oj sin <f>) 



Z s \ — Z\\ 



(15) 



(16) 



where 2h s e is the symmetrical effective length of 
either antenna; 2h a e is the antisymmetrical effective 
length of either antenna ; Z sl is the self-impedance of 
either antenna in the presence of the other, and 
Z 12 is the mutual impedance. Equations (15) and 
(16) follow directly from the equivalent circuit of a 
receiving antenna, for the unloaded condition. 

The no-load currents I\ E (p) and I 2E {o) va^j be 
obtained from (5) and (6) using (15) and (16). 
These are 



I ie (0) = -±ttU + 



cos 



(?*>•)' 



Z s \-\-/ji 2 



+4 



sin 



in ( — sm <i> J 

Zsi Z\ 2 



(17) 
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and 



/ 2 *(0) = 



-4ttUJ 



cos ( — sin <£ J 



Zt s \-\-Zii2 



• h a e 



—3' 



sm ( — sin $ 



'06^ 



A?l — ^12 



(18) 



Second order mutual- and self-impedances are 
available in the literature [King, 1956, ch 3]. If 
the symmetrical and antisymmetrical effective 
lengths can be found to the same accuracy, second 
order no-load currents can be obtained, using (17) 
and (18). It turns out to be much simpler, from the 
computational point of view, to determine the sym- 
metrical and antisymmetrical effective lengths to 
second order accuracy than to evaluate numerically 
the solutions of the integral eqs (7) and (8) to 
obtain second order short-circuit currents. A second 
order solution for the scattering error in direction 
finders is sought in this paper. 

The effective lengths 2h s e and 2h a e are like 2h e for 
an isolated dipole except that the expansion param- 
eters are \{/ s and yp a instead of \p Kl . These differences 
should be expected from an examination of the 
integral eqs (7) and (<S). The various expansion 
parameters are defined as follows when f$h = 7r/2: 



+*= 



c«(i>o)+c b Q,o) 



[King, 1956, p. 268, eqs (12) und (14a)], 
[King, 1956, p. 273, eq (4a)], 



fei = 



a 



(H 



(19) 



(20) 



(21) 



[King, 1956, p. 102, eqs (27) and (29a)]. 
Here 

C tt Q, o)«2 sinh- 1 (J)-Cl w-j Si t, (22) 
C\ Q» o)= Ci fip- Ci pq-j Si fip+j Si fa, (23) 

(24) 



/' 



-V(i)+"4 



-V(i)*+ S! -? 

CiX=0.5772+ln Z-Ci X. 



(25) 



(26) 



An important relation is 

t K1 = n-k 



where 



;> 



-*<?> 



(27) 
(28) 



and £=1.445 when 7 < 12 < 20 for ph=^ [King, 1956, 
p. 105, fig. 20.7]. It is reasonable to suppose that 



and 



xf / a =n a —k. 



[King, 1956, p. 270, fig. 3.3]. 
From (27) and (29), 

and from (27) and (30), 



(29) 
(30) 

(3D 
(32) 



12 s and 12 a may now be computed from (31) and 
(32) using (19), (20), (21), and (28) when f3h=Tr/2. 
The values of \h s e \/\ and \h%\/\ are then read directly 
from figure 2, a curve relating \h e \/X to 12 when 
/3A = 7t/2. It should be mentioned that in general 



h e =K'+jh' e . 



(33) 



Reference to the literature [King, 1956, p. 488, fig. 
9.2a] shows that h' e ~0 for all values of 12 when 
fih = w/2. It is a legitimate procedure, therefore, to 
replace h s e and h a e by \h s e \ and \h a e \, respectively, in (17) 
and (IS) provided the theory is restricted to the case 
ph=T/2. 




Figure 2. Effective length \~\ for a half -wave dipole 
as a function of Q (j3h = 7r/2) 
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In table 1 are listed second order values of Z 12 and 
Z s i as functions of b/\ when fih=ir/2 and 12=10 [King, 
1956, p. 305, tables 8.10 and 8.11]. Table 2 pro- 
vides computed values of Y s = l/Z s and Y a =l/Z a for 
fih=ir/2 and 0=10. This table facilitates error com- 
putations when <J>F^7r/2. Table 3 furnishes computed 
values of \p s , \p a , 12 s , and Q a and the values of \h s e \/\ 
and |A?|/\, as read from figure 2 when 12=10 and 
f3h=w/2. These data permit evaluation of Ii E (o) and 
I 2E (o) from (17) and (18), as functions of <i>, to second 
order accuracy. 



Table 1. Mutual and self-impedance for coupled antennas 



tih - 



Q=2Zn(^)=10, 



(All values are in ohms) 



b 
X 


Zn 


Z,i 


0.1 
.2 
.3 

.4 
.5 

.6 
.7 
.8 
.9 
1.0 


75.1- j 7.8 

55. 2-7 33. 7 

29. 9-7 46. 3 

4.9-/43.2 

-15. 2-> 30. 4 

-25. 2-7 12. 5 
-27.0+7 4.9 
-19.5+7 17.1 
-7. 2+7 21. 6 
5. 1+7 18. 2 


82. 0+7 37. 7 
81. 4+j 38. 7 
82. 0+j 39. 9 

84. 6+j 41. 7 

85. 8+7 42. 7 

86. 2+7 42. 9 
86. 1+j 42. 9 
86. 2+j 42. 8 
86. 3+; 42. 8 
86. 4+j 42. 9 



Table 2. Symmetrical and antisymmetrical admittances for 
coupled antennas 

0fc=|, O=10 



6 
X 


-1 


-1 


0.1 
.2 
.3 
.4 
.5 


0.6143-j 0.1169 
.7311-j 0.0268 
.8909+j 0.0510 
1.112 +j 0.0187 
1. 375 - j 0. 2395 


0. 3258— i 2.149 
.4420—j 1.221 
. 5136— j 0. 8498 
. 5878— j .6261 
. 6498— j .4703 


.6 

.7 
.8 
.9 
1.0 


1. 326 -j 0. 6720 

1.023 -j 0.8273 

0. 8299-j 0. 7453 

. 7603-j 0. 6190 

. 7559— j 0. 5048 


. 7160— j .3529 
. 7946— j .2670 
. 8933— j .2172 
1.017 -j .2307 
1. 126 -j . 3421 



*Multiply all numbers by 10~ 2 



Table 3. Computed values of the expansion parameters and 
the symmetrical and antisymmetrical effective lengths for 

0h=5and 12=10 



b 
X 


4" 


0a 


0« 


Q« 


[A'J/* 


I*:Ia 


0.1 


10.55 


6.78 


12.00 


8.23 


0. 1884 


0. 2077 


.2 


9.18 


8.12 


10. 63 


9.57 


.1944 


.2000 


.3 


8.41 


8.81 


9.86 


10.26 


.1984 


.1964 


.4 


8.03 


9.12 


9.48 


10.57 


.2005 


.1948 


.5 


7.94 


9.17 


9.39 


10.62 


.2019 


.1944 


.6 


8.09 


9.02 


9.54 


10.47 


.2002 


. 1955 


. / 


8.36 


8.77 


9.81 


10.22 


.1988 


.1965 


.8 


8.62 


8.52 


10.10 


9.97 


.1975 


.1979 


.9 


8.81 


8.32 


10.25 


9.77 


.1964 


.1990 


1.0 


8.87 


8.25 


10.32 


9.70 


.1961 


. 1995 



3. Circuital Relations in Coupled Receiving 
Antennas 

Given the short-circuit currents I\ E (o) and I 2E (o), 
the problem is that of finding the voltage V 1L and 
V 2L that is developed across identical load impe- 
dances Z L connected across the terminals in antennas 
1 and 2, respectively. Writing V 1L =\Vi L \/cl> lL and 
V 2I = \V 2L \ Afez,, the angle sought is <£=<£il— </> 2 l. The 
angular error in the interferometer due to antenna 
coupling or scattering is then 



5=(j)—l3b sin $. 



(34) 



The terminal currents in a two-element transmitting 
array immersed in an electric field may be obtained 
by superposition, as indicated in figure 3. Clearly 



and 



/iL(0)=/i*(0) + / ir (0) 

/ 2i (0) = / 2E (0) + / 2V (0) 



(35) 
(36) 



Here Ii E (0) and I 2E (0) are the terminal currents 
when the antennas are excited only by the incident 
electric field, and the impedanceless generator volt- 
ages Vi and V 2 are zero; Ii V (0) and T 2V (0) are the 
terminal currents when the incident electric fields E x 
and E 2 are zero and the antennas are center-driven 
by the voltages Vi and V 2 . 

The fields, currents and voltages shown in figure 3 
are separated into symmetrical and antisymmetrical 
parts as follows: 



E l =E s +E a 
E 2 =E S -E a 
/ li (0)=/i(0)+/£(0) 

/ 2 l(0) = /£(0)-/£(0) 

V 2 =V s -V a 



(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



This separation permits the symmetrical problem 
(superscript s) portrayed by figure 4 to be discussed 
independently of the antisymmetrical problem (su- 
perscript a) portrayed by figure 5. 

Symmetrical Problem 

Refer to figure 4. The following relations are true : 
/J(0j=7i(0)+Jf(0), (43) 

V s 



v s =-pmz l , 

Z s =Z s \ J rZ\ 2 . 



(44) 

(45) 
(46) 
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z = h 



E *t 



• -0 (W) © 

-Yi 2L (o) Y 1 ^' ' 



I 2E (0) 



I 1E (0) 



+ ® , © 



-Ti 2V <o) -yi lv (o) 



z =-h 



Figure 3. The principle of superposition applied to obtain the terminal currents in a 
two-element transmitting array immersed in an electric field. 



-M 



© © 

Vr s L (o) Y i$ l (0) 



ijioi 



l|(0) 



ri»») 



ij») 



Figure 4. 7%e symmetrical problem. 



-°t 






E a t 



I°<0) 



I°<0) 



+Tl°(0) -flS(O) 



Figure 5. 7%e antisymmetrical problem. 
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Hence 



«»='m{zM£z} < 47 > 



Equation (47) establishes the relation between the 
symmetrical part of the load current and the sym- 
metrical part of the short-circuit current. 

Antisymmetrical Problem 

Refer to figure 5. The relations (43) to (45) apply 
with superscript a written for s. The antisymmet- 
rical impedance is required. It is 



It follows that 



Zj a — Zj s \ Zj\ 2 . 



I a M 



=/J(0){ 



Z s \ — Zi< 



Z s \ — Zi 2 -rAc 



(48) 



(49) 



Equation (49) establishes the relation between the 
antisymmetrical part of the load current, and the 
antisymmetrical part of the short-circuit current. 

Figure 6 illustrates the two-element receiving an- 
tenna array with identical impedances Z L connected 
to the load terminals of each antenna. The load 
currents in antennas 1 and 2 are given by (39) and 
(40), respectively. The voltage drop across Z L in 
antenna 1 is then 

V 1L =Z L {Il(0)+I a L (0)}, (50) 

and the voltage drop across Z L in antenna 2 is 

V 2L =Z L {Il(0)-im}- (51) 

ANT I 



ANT 2 



I 2L (0) = lf(0) - I°(0) 



N* 



X |L (0) = I S L (0) + I°(0) 
V |L = I 1L (0)Z L 



Figure 6. Two-element receiving antenna array. 
Impedances of value Zl are connected to the load terminals. 



/£(0) and 7£(0) are expressed in terms of 7^(0) 
and 7J(0) by (47) and (49). Also, from (5) and 
(6), 7i(0) = (7«(0) + / 2if (0))/2 and I%(0) = (I 1B (0) 
— 7 2 #(0))/2. It is now a simple matter to obtain 
V 1L and V 2L , as given by (50) and (51), in terms of 
7i#(0) and 7 2£ -(0). The final expressions are 



V lL = 



(Z sl +Z l2 +Z L ) (Z sl -Z l2 +Z L ) 
{Iie(0)[Z 2 s1 +Z s1 Z l -Z1 2 ]+I 2E (0)Z 12 Z l } 



(52) 



and 



v 2L = 



(Z„+Z n +Z L ) (Z sl -Z 12 +Z L ) 
{I 2E (0)IZ 2 s1 +Z s1 Z l -ZU+Le(0)Z 12 Z l } 



(53) 



These equations provide a way of computing the 
voltages across the identical load impedances in 
terms of the short-circuit currents 7^(0) and 7 2 #(0). 
One check on the validity of either of the above 
expressions is to let b^> oo so that Z 12 ^0. Then (52), 
for example, becomes 



V[ L = 



r 1E (o)Z' sl z L 



AlT^L 



(54) 



Here the primes indicate that the various parameters 
change value when one of the antennas (in this case 
antenna 2) is moved to infinity. Equation (54) is 
recognized immediately as the correct expression for 
the voltage across the load impedance Z L of an 
isolated center-loaded receiving antenna. 
Let the following definitions be made : 



(Z s i-\-Zi 2 -\-Z L )(Z s \— Zi 2 -\-Z L ) 
Z 2 8l +Z sl Z L -Z 2 12 =Aet* 
Z 12 Z L =Be* 
I 2E (0)=CI 1M (0)e" 

With this notation (52) and (53) become 

V^TlMiAe^+BCeW+v] 
V 2L ^Yh{0){ACe j ^ y) l +Be#}- 
Using these expressions and the formula 

/ x-y \ 



(55) 



tan 1 x— tan 1 y=t&n~ 



(56) 

(57) 

(58) 



it is readily deduced that the angle <j> indicated by 
the interferometer is 



<£=tan 



'{ 



(B*-A*)Csm y+AB(l-C 2 ) sin (a- 
(B 2 +A 2 )C cos y+ABil + C 2 ) cos (, 






(59) 
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4. Numerical Analysis 

Figure 7 exhibits a second order curve for 5 as 
a function of fib for an isolated interferometer when 
12 = 10, ph=TJ2, $ = tt/2 and Z L =100+j 0.0 ohms. 
To facilitate computations for other values of <$> and 
Z L , when 12 = 10 and f3h = ir/2, tables 2 and 3 are pro- 
vided. Table 4, giving values of A, a, B and /3 as a 
function of b/\ is useful when constructing curves of 5 
against fib with $ as parameter. This table applies 
only when 12=10, f3h = ir/2 and Z L =100+j 0.0 ohms. 



50 
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O 50 IOO 150 200 250 300 350 400 

ANTENNA SPACING 0b , DEGREES 

Figure 7. Antenna coupling error in a two-antenna inter- 
ferometer. 

fi=10, 0h = ir/2, 3> = 7r/2, and Zl=100+j 0.0 ohms. 
Table 4. A, a, B and versus - 

A 

(These data apply only when ph = irl2, fi=2/w(2/i/a) = 10and Zl=100+jO.O ohms) 



6 
X 


0b 


A* 


a 


B* 


$ 


0.1 


0. 6283 


L.366 


0. 9520 


0. 7551 


-0. 1037 


.2 


1.257 


1.794 


. 8853 


. 6467 


-.547'.) 


.3 


1. 885 


1.974 


. 7396 


.5512 


-.9975 


.4 


2.513 


1. 956 


. 0371 


.4348 


-1.458 


.5 


3.142 


1.826 


.6245 


. 3399 


-2.034 


.6 


3.770 


1.757 


.6783 


.2903 


-2.696 


.7 


4.398 


1.801 


. 7252 


.2744 


2. 962 


.8 


5.027 


1.875 


.7174 


. 2593 


2.422 


.9 


5.655 


1. 893 


. 6849 


.2277 


1. 893 


1.0 


6.283 


1.867 


.6687 


. 1890 


1.300 



5. Concluding Remarks 

It is of interest to observe that when 3>=7r/2, 5=0 
for pb = nw, n = 0, 1, 2, . . . Evidently, when 3>^7r/2 a 
spacing fib exists such that 5 = 0. 

The error curve shown in figure 7 applies to base- 
loaded vertical antennas over a perfectly conducting 
earth when the load resistors are halved. If all im- 
pedances are divided by 2, a, f3, 7, and C remain 
unchanged, but A and B arc one-fourth their free 
space values. Accordingly <j> (and therefore 5) are 
the same as for the isolated array. These statements 
are readily verified using (55) and (59). 

A precise theory has been presented for the scatter- 
ing error in a two-antenna direction finder. It 
applies with equal accuracy to radio interferometers 
employing antennas of length other than A = X/4 
(but with slight increase in complexity). Further- 
more, the theory is readily generalized to include a 
wave tilt with respect to the antenna axes. 



♦Multiply A and B values by 10 4 . 



The writer is indebted to R. W. P. King of Har- 
vard University, and R. H. Duncan of New Mexico 
State University, for reading this paper. He wishes 
to thank Mrs. Dixie 8. Baca for doing the typing. 
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